Carrier-envelope phase dependence in single-cycle laser pulse propagation with the 

inclusion of counter-rotating terms 
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We focus on the propagation properties of a single-cycle laser pulse through a two-level medium 
by numerically solving the full-wave Maxwell-Bloch equations. The counter-rotating terms in the 
spontaneous emission damping are included such that the equations of motion are slightly different 
from the conventional Bloch equations. The counter-rotating terms can considerably suppress the 
broadening of the pulse envelope and the decrease of the group velocity rooted from dispersion. 
Furthermore, for incident single-cycle pulses with envelope area 4-7T, the time-delay of the generated 
soliton pulse from the main pulse depends crucially on the carrier-envelope phase of the incident 
pulse. This can be utilized to determine the carrier-envelope phase of the single-cycle laser pulse. 

PACS numbers: 42.65.Re,42.50.Md,42.50.Ct 
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I. INTRODUCTION 



The modern technological progress in ultrafast optics 
makes it possible to produce few-cycle laser pulses [lJ-Q ■ 
Recently, a single-cycle pulse with a duration of 4.3 fem- 
toseconds has been generated experimentally Fur- 
thermore, great effort for the generation of extremely 
short pulses via few-cycle laser pulses has been made 
[EHI> particularly, single-cycle gap solitons Q and unipo- 
lar half-cycle optical pulses [U , respectively, generated in 
the dense media with a subwavelength structure. If the 
pulse duration approaches the optical cycle, the strong- 
field-matter interaction enters into the extreme nonlin- 
ear optics [9(, and the standard approximations of the 
slowly varying envelope approximation (SVEA), and the 
rotating- wave approximation (RWA) are invalid [l(J [H| ■ 
When the Rabi frequency of the few-cycle laser pulse 
becomes comparable to the light frequency, the electric 
field time-derivative effects will lead to carrier-wave Rabi 
flopping (CWRF) [l2|], which was observed experimen- 
tally in the semiconductor GaAs sample [13 . In this ex- 
treme pumping regime, the simple two-level system can 
still serve as a reference point d, Oj], [H| . 

For the few-cycle laser pulses, the absolute carrier- 
envelope phase (CEP) strongly affects the temporal vari- 
ation of the electric field. These effects give rise to many 
CEP dependent dynamics, such as hi gh- h armonic gener- 
ation [343 i optical field ionization M (20| , atomic co- 
herence and population transfer [2l|, |22| . etc. The CEP 
dependent strong interactions also provide routines to de- 
termine the CEP of few-cycle ultrashort laser pulses. In 
particular, the strong-field photoionization provides very 
efficient tools to measure the CEP of powerful few-cycle 
femtosecond laser pulses for the first time [23[. Another 
promising approach to determine the CEP is introduced 



on the detection of the THz emission by down-conversion 
from the few-cycle strong laser pulse [24|. Recently, the 
angular distribution of the photons emitted by an ultra- 
relativistic accelerated electron also provides a direct way 
of determining the carrier-envelope phase of the driving 
laser field [251 ]. However, all these measurements of CEP 
are based on light amplification in strong-field regime. 

Therefore, it is very meaningful to explore routines for 
determining the CEP of few-cycle laser pulse at relative 
lower intensities without light amplification. The nonper- 
turbative resonant extreme nonlinear optics effects would 
be good candidates for measuring the CEP of few-cycle 
laser pulses with moderate intensities [3, [l5[. However, 
the period of these CEP-dependent effects is tt due to the 
inversion symmetry of light-matter interaction in two- 
level systems. Thus, the sign of the few-cycle laser pulse 
still cannot be determined. In order to remove the tt- 
shift phase ambiguity, the violation of inversion symme- 
try should be considered [26| . In the presence of an elec- 
trical bias, the phase-dependent signal of ultrafast optical 
rectification in a direct-gap semiconductor film implies a 
possible technique to extract the CEP 
the inversion-asymmetry media, such as 1 
[28| and the asymmetric quantum well [2£ 



27|. Moreover, 



solar molecules 
j , could also be 
utilized to determine the CEP of few-cycle laser pulses. 
In this paper, we introduce the counter-rotating terms 
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FIG. 1: (color online) Schematic of the light-matter interac- 
tion scheme. The green curve illustrates the incident few-cycle 
laser pulse. The blue curve illustrates the transmitted laser 
field. Symbol " —> " denotes the propagation direction axis. 
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(CRT) in the spontaneous emission damping, and inves- 
tigate the influence of CRT on the propagation dynamics 
of nonamplihed single-cycle laser pulses in two-level me- 
dia. The CRT should be considered for such ultrashort 
pulses interacting with the medium with strong relax- 
ation processes, because the CRT can notably suppress 
the broadening of the pulse envelope and the decrease 
of the group velocity arising from dispersion. Further- 
more, when the incident single-cycle pulse with envelope 
area Q = 4ir propagates throngh the two-level medium, 
it splits into two pulses. The stronger main pulse moves 
faster than the weaker generated soliton pulse, and the 
pulse time-delay between them shows a pronounced CEP 
dependence. Therefore, in the presence of a static elec- 
tric field, we present a simpler approach for measuring 
the CEP of the few-cycle laser pulses, by detecting the 
time-delay of the generated soliton pulse. 



II. APPROACH 

A. Maxwell Equations 

We consider the propagation of a few-cycle laser pulse 
in a resonant two-level medium along the z axis, as shown 
in Fig. [T] The pulse initially moves in the free-space re- 
gion, then it penetrates the medium on an input interface 
at z — and propagates through the medium, and finally, 
it exits again into the free space through the output in- 
terface at z = L. With the constitute relation for the 
electric displacement for the linear polarization along the 
x axis, D x — eoE x + P x , the full- wave Maxwell equations 
for the medium take the form: 



dHy 

dt 
dE x 
dt 



1 dE x 

Ho dz ' 

1 dH y 1 dP x 



e dz 



eo 



dt 



(la) 
(lb) 



where E x and H y are the electric and magnetic fields, 
respectively. po and eo are the magnetic permeability 
and the electric permittivity in the vacuum, respectively. 
The macroscopic nonlinear polarization P x = —Nd%2U is 
connected with the off-diagonal density matrix element 
P12 = h(u + iv) and the population inversion w — P22 — 
pn, which are determined by the Bloch equations below. 



B. Master Equation 

The Hamiltonian of the two-level system we considered 
can be described by [30j : 

H = ^hw k a\a k + hoj S z + hn(t)(S + + S~) 

k 

+ i^2(g k -d 21 ){al(S + + S-)-R.c.}, (2) 



where loq is the transition frequency, and c?2i is the elec- 
tric dipole moment of the transition between the upper 
state |2) and the lower state |1). at (<Zfe) is the creation 
(annihilation) operator for photons with momentum hk 

and energy foujk, while gf. = ^/ 27r ^" fc e\ describes the 
vacuum-atom coupling and e\ represents the unit polar- 
ization vector with A e {1,2}. S+ = |2)(1| (S~ = |1}(2|) 
is the dipole raising (lowering) operator of the two-level 
system, S z — (|2)(2| — |l)(l|)/2 is the inversion operator. 
f2(t) — dvzEx/K is the Rabi frequency of the incident 
laser field. 

In the usual Born-Markov and mean-field approxima- 
tion, but without the rotating-wave approximation, the 
master equation of the system is determined by 



P (t) + i[u s z + n(t)(s + + s-), P ] 

= - 1 [S+,(S+ + S-)p] +H.c, 



(3) 



where an overdot denotes differentiation with respect to 
time. Here, [S + , S + p(t)] and its hermitian conjugate 
term represent the counter-rotating terms (CRT) for the 
spontaneous emission damping, which are neglected un- 
der the rotating-wave approximation when the duration 
of the laser field pulse r p is much larger than ojq 1 . How- 
ever, for the few-cycle pulses, even the single-cycle or 
sub-cycle pulse, the CRT become indispensable and can- 
not be neglected. In the following, we will investigate the 
effects of the CRT on the propagation dynamics of the 
single-cycle laser pulse in the two-level medium. 



C. Bloch Equations 

Based on the master equation ([3]), including the CRT 
in the spontaneous emission damping, the Bloch equa- 
tions with CRT can be easily derived as follows: 



u = ujov, 

v = —ujqu + 2Q(t)w — 272V, 
w= -2f2(t)t!-7i(io + l), 



(4a) 
(4b) 
(4c) 



where 71 and 72 are the spontaneous decay rates of the 
population and polarization, respectively. The Bloch 
equations with CRT [Eqs. (UJ)] are slightly different 
from the conventional Bloch equations (see for instance 
Refs. HE3): 



U = OJ V - 72U, 

V = —LdQU + 2Q(t)w — 72V, 

w = -2tl(t)v - 71 (w + 1), 



(5a) 
(5b) 
(5c) 



in which the relaxation constants 71 and 72 are added 
phenomenologically. 



D. Numerical Method 

The propagation properties of the few-cycle laser pulse 
in the two-level medium can be modeled by the full-wave 
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the vector potential at z 
A x (t) = A Q sech[1.76(i- 



= as in Refs. [32], |3 

to)/r p ] sin [u p (t - t ) + 4>], (6) 
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FIG. 2: (color online) (a) The Rabi frequency of the incident 
single-cycle pulse with envelope area = 2n. (b) and (c) 
are the time-dependent electric fields and the corresponding 
population inversions at the distance z = 90 fim, respectively. 
The length of the two-level medium is choosen as L — 110 /um. 
The blue solid curves are for the case of Maxwell-Bloch equa- 
tions with CRT, while the red dashed curves are for the case 
with the conventional Maxwell-Bloch equations. The black 
lines represent the pulse envelope. 



Maxwell-Bloch equations beyond the SVEA and RWA, 
which can be solved by the iterative predictor-corrector 
finite-difference time-domain discretization scheme (TTI . 
l3lj . For such an extremely short laser pulse, we define 
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FIG. 3: The corresponding spectra of the electric field in 
Fig. HJb). The black dashed-dotted curve is the spectrum 
of the incident laser pulse. The blue curve depicts the case 
with CRT, and the red dashed curve the case without CRT. 



where Aq is the peak amplitude of the vector potential, 
u),p is the photon energy, and <f> being the CEP. r p is the 
full width at half maximum (FWHM) of the short pulse 
and to is the delay. The electric field can be obtained from 
E x = —dA x (t)/dt. In what follows, we assume that the 
two-level medium is initialized in the ground state with 
u — v = and w = — 1. The material parameters are 
chosen as in Ref. [5j: ojq = 2.3 fs _1 (A = 830 nm), d\2 — 
2 x 10~ 29 Asm, 7f 1 = 1 ps, 7^ = 0.5 ps, the density 
N = 4.4 x 10 20 cm" 3 . The incident pulse has a FWHM 
in single optical cycle t v — 2.8 fs and the photon energy 
u! p = uiq. The Rabi frequency Qq = —A$u} p d/h = 1 fs -1 
corresponds to the electric field of E x = 
an intensity of I = 6.6 x 10 12 W/cm 2 



5 x 10 9 V/m or 
and the incident 



pulse area is defined as = /"L Q(t)dt. 



III. RESULTS AND DISCUSSION 

Now, we focus on the effects of CRT on the propagation 
dynamics of single-cycle laser pulses in two-level medium 
by comparing the numerical results from the Maxwell- 
Bloch equations with CRT [Eqs. (p} and (j4|)] and without 
CRT [Eqs. (p} and ©]. We use an incident single-cycle 
pulse with envelope area = 27r for these simulations 
with the medium zone length: z = 110 /im. 

According to the standard area theorem, the pulse with 
area = 2-7T can transparently propagate through the 
two-level medium without suffering significant lossness - 
the so-called self-induced transparency (SIT) [34|. How- 
ever, when the laser pulse envelope contains only few op- 
tical cycles, the standard area theorem breaks down be- 
cause of the occurrence of CWRF [12| . From our numer- 
ical results, for the short propagation distance, the usual 
SIT regime is essentially recovered. However, at a further 
distance, the established conditions for SIT are destroyed 
due to the extreme nonlinear optical effects. Fig. GJb) 
and Fig. [2jc) present the normalized electric-field pulses 
and the corresponding population inversions at the dis- 
tance z = 90 ^tm for different approaches, namely, the 
blue solid curves depict the case obtained from Maxwell- 
Bloch equations with CRT, while the red dashed curves 
are for the conventional approach without CRT. Com- 
pared with the incident single-cycle 27r pulse in Fig. [2][a) , 
the electric-field pulses for both two cases become clearly 
broadened induced by the dispersion, and suffer the de- 
crease in pulse amplitude. Accordingly, the population 
differences for both cases undergo an incomplete Rabi 
flopping with the CWRF [Fig.^c)]. 

However, there are notably different features between 
these two approaches. The electric-field pulse from the 
approach with CRT [blue solid curves in Fig. Hfb)] is ev- 
idently narrower than that in the case without CRT [red 
dashed curves in Fig. |2Jb)]. This can be easily seen from 
the corresponding spectra shown in Fig. [3] The spectrum 
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FIG. 4: The transmitted electric field of the incident pulses 
with envelope area <d — <in for different initial CEPs: (a) 
and (b) n/2. The length of two-level medium is L — 80 fim. 



of the case with CRT is obviously more broadened than 
that in the case without CRT, although both of them 
become narrower than the incident spectrum. The en- 
velope peak from the approach with CRT is relatively 
larger than that in the case without CRT [Fig. EJJb)], 
hence, the former case lead to more population inversion 
at the leading edge of the electric-field pulse [Fig. [2jc)]. 
Moreover, there is a notably time delay of the electric- 
field pulses and the corresponding population inversions 
between the two approaches [Fig.[5Jb) and (c)]. It means 
that the group velocity of the propagating pulse from the 
conventional approach without CRT is obviously smaller 
than that in the case with CRT. This difference in the 
group velocity is rooted from the different influence of 
dispersion effects for these two approaches. Compar- 
ing Eqs. Q with Eqs. ([5]), there is no damping of the 
real part of polarization u in the Bloch equations with 
CRT, which indicates that the dispersion does not suffer 
lossness. That is to say, the presence of CRT evidently 
suppresses the strong dispersion effects, which lead to 
the broadening of pulse envelope and the decrease of the 
group velocity. 

In addition, we also find that the influence of CRT on 
the propagation dynamics of the single-cycle laser pulses 
is significantly enhanced with the increase of the sponta- 
neous decay rates. Therefore, the CRT is important and 
indispensable for the study of the propagation properties 
of few-cycle laser pulses in the medium with strong relax- 
ation processes. In the following discussion, we will use 
our established full-wave Maxwell-Bloch equations with 
CRT [Eqs. ([T]) and @] to explore an approach for deter- 
mining the CEP of the single-cycle laser pulse. 

In what follows, we simulate the incident single-cycle 
pulses with larger envelope area, i.e. = 47r, prop- 
agating through the two-level medium with a length 
L = 80 /im. During the course of pulse propagation, 
the medium absorbs and emits photons and redistributes 



energy in the pulse. The propagating pulses are altered 
in shape until it reaches a stable status after some prop- 
agation distance by splitting into two pulses, the strong 
main pulse and the SIT soliton pulse. However, the for- 
mer moves faster than the latter, which is why the gener- 
ated SIT soliton pulse breaks up from the main pulse. We 
show the transmitted pulses of the incident single-cycle 
pulse with pulse envelope area = 47r for different CEP 
<f> = and <j> — ir/2 in Fig. [4] It can be seen that both of 
the transmitted pulses split into two pulses. There is a 
time delay between the main pulse and the soliton pulse 
defined as t((f>). The time-delay for the incident pulse 
with CEP cf> = ir/2 [t((f> — tt/2)] is evidently larger than 
that in the case with CEP <f> = \t{cf> = 0)]. It demon- 
strates that the pulse time delay t(4>) is sensitive to the 
initial CEP of the incident pulse. 

For simplicity, we define the relative pulse time delay 
At = t(cf>) — t(cf> = 0) to indicate the CEP dependence. 
We present the relative pulse time delay At as a function 
of the initial CEP of the incident pulse in Fig. [S] with blue 
circles. It is found that the relative pulse time delay At 
is related to the CEP of the incident pulse with a nearly 
cosinelike dependence. However, the time delay t((f> = tt) 
is exactly the same as t(<fi — 0), and hence, the period of 
the CEP-dependent pulse time delay is only n because of 
the inversion symmetry of light-matter interaction. This 
means that we cannot distinguish the incident pulse from 
the initial CEP <j> -)• cj) + n. 

In order to remove the 7r-shift phase ambiguity, we 
add a static electric field to break the inversion sym- 
metry of the light- matter interaction [35] ]. As a result, 
the Rabi frequency terms in Bloch equations [Eqs. ((H)] 
should change as ft(t) — > il(t) + f, where / describes the 
strength of the static electric field. The presence of the 
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FIG. 5: The relative pulse delay between the transmitted 
soliton pulses as a function of the CEP of the incident pulses 
for different strengths of the static electric field. The solid 
lines are a guide for the eyes. 
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static electric field gives rise to the enhancement of the 
CEP-dependent variation in the peak electric strength 
of the single-cycle pulse, which will enhance the CEP 
dependence of the dynamics effects. The relative pulse 
time delay At of the transmitted soliton pulses as a func- 
tion of the initial CEP of the incident single-cycle pulses 
for different static electric fields is presented in Fig. [5] 
Compared with the blue circles of / = 0, the influence 
of the static electric field on the relative pulse time delay 
is significant. Let us take / = 2% Q , for example green 
squares in Fig. [5j then the relative time delay At ^ 
at = 7T, i.e., t((f> = tt) is quite different from t{<j) = 0) 
in the presence of the static electric field. The variation 
with the CEP of the incident pusle also becomes much 
stronger. The period of the CEP-dependent pulse time 
delay becomes 2n because the inversion symmetry is bro- 
ken assisted by the static electric field. Moreover, with 
the increase of the static electric field, such as / = 3% flo 
and / = 4% f2 , the dependence of the relative pulse time 
delay on the initial CEP is further enhanced [red diamond 
and black circles in Fig. [5] . 

As a result, in the presence of the static electric field, 
if the relative time delay of the generated soliton pulses 
is calibrated, this effect suggests an approach for deter- 
mining the CEP of the incident single-cycle laser pulses 
in both sign and amplitude. In addition, it should be 
pointed out that the pulse time delay might be much eas- 
ier to detect compared with other features of the soliton 
pulse, such as the intensity and pulse duration [28[. Fi- 
nally, in our discussion, the static electric filed strengths, 
which are a few percentages of the single cycle laser pulse 
strength, exceed a few MV/cm. In order to achieve this 



kind of strength of the static electric field in an experi- 
ment jWe may proceed with a special case as suggested in 
Ref. [35[ where an additional electric field with a much 
lower frequency (such as CO2 laser field, terahertz pulses 
or a midinfrarcd optical parameter amplifier pulse) is 
used instead of the static electric field. The ultra-short 
dynamics can prevent the system from being destroyed 
or ionized. 



IV. SUMMARY 

In summary we investigated the propagation proper- 
ties of single-cycle laser pulses in a two-level medium in- 
cluding the counter-rotating terms in the spontaneous 
emission damping. We found that the counter-rotating 
term can efficiently suppress the broadening of the pulse 
envelope and the decrease of the group velocity. Thus, 
the counter-rotating term is important and indispensable 
for the study of the propagation dynamics of few-cycle 
laser pulses, even for single-cycle and sub-cycle pulses. 
Furthermore, we explored the CEP-dependence of the 
generated soliton pulse from the single-cycle laser pulse 
propagating through the two-level medium. The time 
delay of generated soliton pulses depends sensitively on 
the CEP of single-cycle incident laser pulse. Hence, the 
presence of the static electric field enhances the CEP- 
dependence of the relative pulse time delay, which have 
an potential application in determining the CEP of the 
incident single-cycle laser pulse. 
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